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We determine conditions under which a spin-liquid Mott insulator |0) defined by a Gutzwiller 
projected BCS state at half-filling is Z2 fractionalized. We construct a trial vison (Z2 vortex) 
state \V) by projecting an hc/2e vortex threading the hole of a cylinder/torus and examine its 
overlap with |0) using analytical and numerical calculations. We find that generically the overlap 
vanishes in the thermodynamic limit so the spin-liquid is Z2 fractionalized. However, for microscopic 
parameters appropriate for high Tc cuprates, we estimate the vison gap E v <C J, consistent with 
recent experimental bounds, due to the proximity to the bipartite symmetric point where E v = 0. 
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Spin liquid insulators doped with holes have long been 
regarded as natural candidates for high temperature su- 
perconductivity since the resonating valence bond (RVB) 
proposal of Anderson Q . There has been a revival of in- 
terest in this idea since such wavefunctions have recently 
been shown to provide a remarkably good description |2| 
of the superconducting (SC) state of the high Tc cuprates 
over a wide doping range. In addition, spin liquid states 
for frustrated magnets are of great interest both theoret- 
ically 3] and experimentally 

It has also long been suspected that RVB insulators 
may exhibit quantum number fractionalization or spin- 
charge separation i.e., support neutral, spin- 
1/2 and charged, spin-0 excitations. Senthil and Fisher 
[lfj have recently emphasized that such fractionalized 
states in d=2, 3 may be described as deconfincd phases of 
Z2 gauge theories, and necessarily imply a novel gapped 
excitation jy, la, : the Z2 vortex or "vison" . The pres- 
ence/absence of visons threading noncontractible loops 
leads to a topological ground state degeneracy on cylin- 
ders or torii They further proposed an inge- 
nious flux-trapping experiment to detect the vison. How- 
ever, a series of experiments to detect visons in highly 
underdoped cuprates have led to negative results, setting 
bounds on the vison gap E v < 150— 180K 12j, an order- 
of-magnitude smaller than the natural scale J ~ 10 3 K. 

Motivated by the above developments we address the 
following issues: (a) How can one think about the 
vison in terms of electronic coordinates? (b) Under 
what conditions is the vison well-defined, leading to Z2- 
fractionalization and topological degeneracy in spin liq- 
uid insulators? (c) How large is the vison gap, the en- 
ergy scale below which fractionalization is apparent? (d) 
What are the implications for cuprates and for other frus- 
trated magnetic materials? 

Vison Wavefunction: For conceptual clarity, consider 
first a square lattice of L 2 sites "wrapped" into a cylin- 
der along the x direction. The spin-liquid ground state 



0) of interest to us is given by the Gutzwiller projec- 
tion V = Y\ r (l — n r ^n T i) of the Af-particle d-wave BCS 
wavefunction at half-filling (N = L 2 ): 

|0) = V\BCS) = 7>[5>(r - r')4 T 4,J W/a |vac>. (1) 

r.r' 

Here <p{p) = L~ 2 £ k exp(zk • p) [A k /(& + v^k + A k)] » 
the internal pair wavefunction with the variational pa- 
rameters |l4j {j, and A which determine £k = e(k) — fi, 
with e(k) = — 2t (cos k x + cos k y ) — At 1 cos k x cos k y and 
Ak = A (cos fez — cos fey) /2. The k's are chosen consis- 
tent with periodic boundary condition (PBC) along x. 

With optimized variational parameters, V\BCS) is an 
energetically good variational ground state for the large- 
U Hubbard model |2, ll3j- However, we focus more 
generally on understanding topological degeneracy and 
fractionalization in the phases described by this class of 
states rather than just the energetically optimal wave- 
function. Toward this end, we define a trial vison state 
\V y ) by Gutzwiller projecting an hc/2e vortex threading 
the cylinder along p by \V y ) = V\ (hc/2e) ), follo wing; 
Senthil and Ivanov [15j and building on earlier ideas [la] • 

To construct the hc/2e vortex we simply modify the 
pair cp(r — r') in the BCS ground state to expfiQ • (r + 
r')/2]iy9 J 4(r—r'), where Q = 2ttx/L twists the phase of the 
condensate by 27r going around the cylinder. However, 
the requirement of a single-valued wavefunction, when r 
(or r') goes around the cylinder, constrains k's to sat- 
isfy antiperiodic boundary conditions (APBC) along x 
in the Fourier transform defining ^(r — r'); hence the 
subscript A. For the fully projected state at half-filling 
one can show that the phase factors drop out; in first 
quantized form, these are independent of the spin con- 
figuration. Then the trial vison state \V y ) simplifies to a 
form identical to eq. JQ) with ip — > tpA- 

Note that a projected hc/e vortex with Q = 4ttx/L 
does not require any change in boundary conditions on 
ip and is therefore identical to the ground state [0). The 
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FIG. 1: (A) "Phase diagram" showing that P\0) is Z 2 
fractionalized in the entire (£', u)-plane except at the bi- 
partite symmetric point (0,0) (see text for details). (1) 
Shaded and unshaded half-planes are related by symmetry 
(t',n) — * (—t,—fi); (2) hatched region for large negative fi 
is the short-range RVB limit; (3) analytical arguments for 
/i = ^ibcs (schematic line) are described in the text; (4) 
crosses indicate some points where we numerically compute 
overlaps; (5) shaded ellipse around (—0.25,0) is the region 
relevant for high Tc SCs. (B) Energy of |0) for the Hub- 
bard model at half-filling for parameters relevant to cuprates 
(U/t=12, t'/t=-0.25 and optimal A=1.25f) versus u/t. The 
optimal value is —0.3 < fi/t < 0.3 within statistical errors. 

manifest Z2 character (two visons is the same as none) of 
projected vortices makes them attractive candidates for 
visons. On a torus |0) corresponds to PBC along both 
x and y, the vison states \V X ) and \V y ) correspond to 
PBC/APBC and \V xy ) to APBC along both directions. 

Phase diagram and Symmetries: A superconduct- 
ing state with a vortex threading the hole of the cylin- 
der/torus carries current and is orthogonal to the ground 
state. However, Gutzwiller projecting the BCS wavefunc- 
tion at half-filling destroys SC order |2| and results in a 
spin-liquid insulator |l7|. The question we would like to 
address is: does any remnance of vorticity survive projec- 
tion? If it does, then (V a \0)={V a | Vg)=0, for afi—x, y, xy, 
the vison is well-defined leading to topological degener- 
acy |l8j . and one is in a Z2 fractionalized phase. If, on the 
other hand, the proposed vison state has nonzero overlap 
with the state |0) in the thermodynamic limit, one is in 
a conventional non-fractionalizcd phase. Our main task 
then is to compute these overlaps and establish the exis- 
tence of disconnected topological sectors as evidence for 
fractionalization. 

To determine the conditions under which the vison sur- 
vives, it is useful to consider the full parameter space for 
projected BCS states at half-filling, which is the (t',fi) 
plane with t = 1 and A held fixed at some non-zero 
value in order to describe a RVB liquid of singlet pairs. 

We can use symmetry arguments to show that not 
all states in this space are distinct upon projection and 
|0(t',/z)) = |0(— t',— fi)). In brief, we can change t' — > 



—t 1 , fj, — > — fi by a global particle-hole transformation 
c ra — > (— l) x+v cl a in the wavefunction, also redefining 
the I vac) since empty sites transform into doubly occu- 
pied ones. However, with exactly one particle per site 
such a transformation interchanges f-spins and j-spins 
at every site for any configuration. Since the BCS wave- 
functions are spin-singlets this leaves the state invariant. 
We thus restrict attention to t' < in Fig.QjA). 

Short-range RVB limit: Note that large neg ative \i 
with [i < — 4(|t| + \t'\) is the bosonic limit |l9j of the 
BCS wavefunction, with y(r — r') decaying exponentially 
in real space. In this limit |0) may be viewed as a short- 
range RVB state. Earlier studies of closely related dimcr 
models H,@ indicate that, for short-range singlet bonds 
there exist four topological sectors on a torus, which may 
be straightforwardly related to sectors with/without vi- 
sons in the x, y directions. Further, two dimer configura- 
tions from different sectors have exponentially small over- 
lap ~ exp(— aL). A state is a superposition of many con- 
figurations, and assuming that individually small over- 
laps do not add up coherently, the overlap of two states 
from different sectors would vanish as L — > 00. 

We have numerically checked (see below) that the over- 
lap vanishes for large L validating the above assumption, 
and implying fractionalization in the short range RVB 
limit. Using the symmetry discussed above, we conclude 
that states in the region with (positive) fi > 4(|t| + \t'\) 
are also fractionalized, even though <p(r — r') is not obvi- 
ously short-ranged here. 

Analytical insights for fi = hbcs'- We next focus on 
states for which fi = ^BCs(i') (see the curve in Fig.[TjA)), 
where /iecs is the value at which the unprojected BCS 
state is at half-filling. The unprojected state may then be 
viewed as the ground state of a d-wave BCS mean field 
(MF) Hamiltonian and it is a coherent superposition of 
number eigenstates, sharply peaked at the correct mean 
density. Gutzwiller projection picks out the N = L 2 
contribution with no double occupancy. 

Define local SU (2) gauge transformations [2(J U, gen- 
erated by T+ = ct T c| 4 , ZT = Cil C ih 77 = 4a C ia ~ X ' 

which mix empty and doubly occupied sites as an SU (2) 
doublet, but act trivially on the projected subspace with 
exactly one particle per site. Gutzwiller projection is then 
equivalent to projection onto the SU(2) singlet subspace, 
and we may write any state V\$) = f u U\Q), where the 

integral is over all group elements Ug — exp(i J2 i Ti ■ 6i). 
The overlap of interest can then be written as 

(V|0) = / {{hc/2e)\U\BCS), (2) 
Ju 

which now reduces the problem to computing overlaps of 
unprojected states. 

For a nonbipartite lattice (e.g., with t' ^ 0), one can- 
not gauge away the off-diagonal term in the d-wave BCS 
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MF Hamiltonian and U\BCS) is a SC state (or vortex 
vacuum) for arbitrary U 15]. Then ((hc/2e)\U\BCS) = 
for arbitrary U and (V|0) vanishes. Thus non- 
bipartiteness is a sufficient condition for fractionalization 
in d-wave states when [i — /ibcs H3- 

On the other hand, for a bipartite system (in our case 
t' = 0), it is well known [i(J that one can gauge away 
the off-diagonal part of the MF Hamiltonian (provided 
it doesn't break time reversal) and transform the BCS 
state into a staggered flux state. Now {{hc/2e)\U\BCS) 
is non-zero for some choice U and one cannot use eq. (j2J 
to argue that (V\0) vanishes. 

Overlaps within Z 2 gauge theories: Before proceed- 
ing to the numerical results, it is useful to understand 
how overlaps like (V\0) are expected to scale with sys- 
tem size L 2 . In a non-fractionalized phase the overlap 
will remain non-zero as L — ► oo. However, in a frac- 
tionalized phase we argue that it should vanish exponen- 
tially for large L. This is easy to see when the matter 
fields are gapped (e.g., the short-range RVB limit) and 
may be integrated out to obtain a 2 + 1 dimensional Z 2 
gauge theory. Deep in the deconfining phase of a pure 
gauge theory, we have evaluated the overlap perturba- 
tively and it scales as (V|0) ~ exp(— L/£). To under- 
stand this result note that the overlap in the deconfined 
phase behaves similar to the spin-spin correlation of the 
dual fsing model at distance L, where the dual spins rep- 
resent the gapped visons. Since the dual Ising model 
is disordered when the gauge theory is deconfining, the 
exponential decay follows. Another way to understand 
this result is that the overlap is nonzero due to tunneling 
of the vison (a "string" of length L) out of the hole of 
the cylinder. Such tunneling events involve intermediate 
states in which the vison pierces the surface of the cylin- 
der. This costs an energy proportional to the vison gap, 
and leads to an overlap exponentially small in L. The sit- 
uation with gapless matter fields is less well studied, even 
though numerically (see below) we find clear evidence for 
exponential decay. 

Numerical results: We finally discuss the results of nu- 
merical evaluation of vison overlaps using the variational 
Monte Carlo method. For computational simplicity j2]| 
we focus on the overlap (V x \V y ). We choose t — 1 and 
fix A = 1.25 for all the calculations reported below. The 
following numerical results are obtained in regimes where 
we have given analytical arguments above, (f) We have 
explicitly checked that the overlaps are invariant under 
(t' , /i) — ► (— t', — /j,). (2) In the regime of short-range RVB 
we find that the system is Z 2 fractionalized; see Fig. 2(A) 
for large negative [i. (3) On the curve [i = /ibcs we find 
that the non-bipartite systems are fractionalized. (4) For 
the special point (t' , /i) = (0, 0) with bipartite symmetry, 
we indeed find that the overlaps are non-zero, indepen- 
dent of L (see Fig. 2(A)), the vison is not well defined 
and the system is not fractionalized. Although consis- 
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FIG. 2: Size dependence of the overlap {V x \V y )(L) for A = 
1.25 with (A) t' = and (B) t' = -1/4. The various curves 
correspond to the /j, value indicated. The lines are fits to the 
form (Vx\Vy)(L) = a(l — tanh((L — which works well 

for the entire regime of parameters we have studied, and from 
which we extract £* . (C) Exponential asymptotic behavior of 
the overlap shown on semi- log scale. 

tent with the analytic arguments, the numerical results 
are important because the arguments were not rigorous 
for (2) and (3), and silent on point (4). 

We next evaluate overlaps at t' = 0.0,-0.25,-0.5 for 
— 1.5 < M ~ 1-5 which, as we discuss below, covers the 
region of interest for possible spin-liquid insulators in the 
vicinity of the high Tc cuprates. We show in Fig. the 
overlap (V x \V y ) as a function of system size L for a range 
of t' and \i. We find that at small L the overlaps are finite 
and then cross over on a scale £* to an asymptotic decay 
cxp(— 2L/£). This behavior can be simply described by 
the functional form (V x \V y ){L) = a(l -tanh((L -£*)/£)) 
which appears to fit the data well. (The factor of two is 
included in the asymptotic decay definition since we are 
looking at overlap of states with two visons in different 
directions). In the regime of parameters where we can 
reliably extract both and £, we find £* ~ 2£. More 
generally, it is hard to access the asymptotic behavior in 
the region of interest, and we can extract crossover scale 
£* more accurately than £. 

We plot the inverse length scale l/£* for t' = 
0, -0.25, -0.5 at different /i in Fig. For t' = 0, we 
see that £* — * oo as fi — > 0, fully consistent with the fi- 
nite overlap independent of L on accessible system sizes 
for (t',fj,) = (0,0) seen in Fig.^A). It also clear that C 
remains finite everywhere away from the origin though it 
may become quite large in its vicinity. 

Since fractionalization manifests itself at length scales 
larger than £*, our numerical results strongly suggest 
that all points in the phase diagram of Fig. QA) are 
fractionalized, except for the origin where £* diverges. 
(t',/i) = (0,0), with its special bipartite symmetry, is a 
non-fractionalized "singular point" in the space of spin- 
liquid insulators that we study. 

Vison gaps: Now the question arises: what regime in 
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FIG. 3: Dependence of l/£* on fi for (A) t' = 0.0, (B) t' = 
—0.25, (C) t' = —0.5. In panel (A), l/£* extrapolates to zero 
on approaching (t , fi) = (0, 0) indicating a non-fractionalized 
"singular point" at the origin. For t' = —0.25, —0.5 (B,C) the 
length scale appears finite everywhere. 

parameter space is relevant for the cuprates? One pos- 
sibility is that there are no fractionalized states in the 
vicinity of the observed phases in real materials. An 
alternative worth exploring in view of the success of 
RVB wavefunctions in understanding the SC state 0, 
is that we take the insulating limit (hole doping x— >0) 
of V\BCS). Optimizing variational parameters /x, A for 
the large-C/ Hubbard model (£7= 12,*' = -1/4), we find 
A~1.25 while -0.3 < fi < 0.3 as shown in Fig.^B). For 
this region (the shaded ellipse in Fig.[TjA)) we conclude 
from Fig. IHlB) that £* > 25 lattice spacings. 

We now convert this to an estimate of the vison gap 
E v . We expect E v = aJ/(C) z < aJ/£,* which van- 
ishes at the "singular point" with bipartite symmetry. 
Here J is the nearest-neighbor superexchange, a = 0(1) 
is a dimensionless constant and the dynamical exponent 
z>\. For a=l and J=1200iV the estimated E v < 50K 
(but this energy could be higher if a is larger). More 
concretely, we emphasize that proximity to the 'bipar- 
tite point' (t' = 0, = 0) can lead to a very small vison 
gap for fractionalized RVB states. This is consistent with 
existing experimental bounds |l2]| and suggests that frac- 
tionalization is apparent only at very low T and does not 
play a role in pseudogap anomalies. 

The ideas and methods developed here should be ex- 
plored further in highly frustrated magnets, with large 
t' or ring-exchange interactions |23| , which may be more 
promising candidates for observing fractionalization. 
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